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ABSTRACT 



Objective MicroRNA plays a vital role in gene expression, and microRNA dysregulation is involved in carcinogenesis. The miR- 
196a-2 polymorphism rsl 1614913 is reportedly associated with cancer susceptibility. This meta-analysis was performed to assess the 
overall association of miR-196a-2 with cancer risk. 

Methods A total of 27 independent case-control studies involving 10,435 cases and 12,075 controls were analyzed for the rsll614913 
polymorphism. 

Results A significant association was found between rsll614913 polymorphism and cancer risk in four genetic models (CT vs. TT, 
OR=1.15, 95%CI=1. 05-1.27; CC vs. TT, OR=1.23, 95%CI=1. 08-1.39; Dominant model, OR=1.17, 95%CI=1.06-1.30; Additive model, 
OR=1.08, 95%CI=1.01-1.14). In the subgroup analysis of different tumor types, the C allele was associated with increased risk of lung, 
breast, and colorectal cancer, but not with liver, gastric, or esophageal cancer. In the subgroup analysis by ethnicity, a significantly 
increased risk of cancer was found among Asians in all genetic models, but no associations were found in the Caucasian subgroup. 
Conclusions The meta-analysis demonstrated that the miR-196a-2 polymorphism is associated with cancer susceptibility, especially 
lung cancer, colorectal cancer, and breast cancer among Asian populations. 
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Introduction 

Cancer, as a major public health problem, is reportedly 
one of the leading causes of death worldwide m . Moreover, 
cancer is a very complex genetic disease, the mechanism 
of which has not been completely elucidated. Studies have 
stiggested that cancer development results from gene- 
environment interactions [2 '. The presence of high-frequency 
low-penetrance susceptibility genes may be important in 
carcinogenesis. 

MicroRNA (miRNA) are 21 to 24 nucleotide-long, single- 
stranded, non-coding RNA. These evolutionarily highly 
conserved miRNA play vital roles in the regulation of 
gene expression through mRNA cleavage or translational 
repression [3 - 41 . To date, 678 human miRNA have been 
characterized; they regulate the expression of approximately 
one-third of human genes' 5,61 . miRNA are possibly involved 
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in many biological processes, including cell differentiation, 
proliferation, and apoptosis p] . Evidences also suggest that 
miRNA may play important roles in tumor development and 
prognosis t8,9) . 

Two types of single nucleotide polymorphisms (SNPs) 
were discovered in miRNA genes based on their location. 
One type is located in the mature regions of the miRNA and 
the other is located in the pre-miRNA regions. The former 
may directly regulate both the binding to target mRNA 
and pre-miRNA maturation, whereas the latter may only 
influence the maturation of pre-miRNA [ia n \ 

The SNPs in pre-miRNA have drawn increasing attention 
because they influence the maturation of miRNA, and they 
play potential roles in tumor development and progression 1121 . 
In 2008, Hu et al. 1131 identified a polymorphism in miR-196a-2 
with a T to C change (rsll614913). Rsll614913 is located in 
the 3' passenger (3p) strand mature sequence of miR-196a-2. 
This functional polymorphism is reportedly associated 
with the susceptibility of various tumors, including lung 
cancer [14 - 151 and breast cancer ,161 , and with lower survival 
rates of non-small cell lung cancer [13] , gliomas [17] , gastric 
cancer [1S1 , gallbladder cancer [w \ head and neck cancer [201 , 
esophageal cancer [211 , and hepatocellular carcinoma |22 ' 231 . 
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Since then, many epidemiologic studies have examined the 
associations between SNPs and the risks of various cancers 
in diverse racial/ethnic populations. However, the results of 
these studies are inconclusive because of their small sample 
size or lack of replication. Therefore, this meta-analysis was 
employed to investigate the association of the miR-196a-2 
rsll614913 polymorphism with cancer risk. 

Materials and Methods 
Literature search 

A meta-analysis was conducted based on all case-control 
studies that examined the association between miR-196a2 
rsll614913 polymorphism and cancer risk. The terms 
rsll614913, polymorphism, and cancer were used to search 
in the PubMed, EMBASE, Cochran library, and Web of 
Knowledge databases. The last search was completed on 
January 31, 2012. No language restrictions were imposed in 
the meta-analysis. Two independent investigators performed 
the search, Zhang H and Su YL. The studies that fulfilled 
the inclusion criteria were selected for evaluation. A manual 
search was also conducted through reviewing the references 
cited by the selected articles. 

Inclusion criteria 

The inclusion criteria for our meta-analysis were as follows: 
(1) a case-control study design; (2) an association reported 
between miR-196a2 rsll614913 polymorphism and cancer 
risk; (3) number of subjects in each genotype available or 
sufficient allele frequencies for estimating odds ratio (OR), 
and 95% confidence interval (CI). 

Data extraction 

Two independent investigators evaluated all the studies 
selected and extracted the data as shown in Table 1. When a 
discrepancy was found, a third investigator was invited for 
discussion until a consensus was reached. The extracted data 
included the name of the first author, year of publication, 
country origin, ethnicity, cancer type, source of controls 
(population- or hospital-based), total number of cases and 
controls, and the number of different genotypes in the 
respective cases and controls. 

Statistical analysis 

The consistency of the data with the Hardy-Weinberg 
equilibrium (HWE) was determined using a % 2 test after 
comparing the observed and expected genotypic frequencies 
in the control populations; a P-value <0.05 was considered 
statistically significant (Table 1). 

The association between miR-196a-2 polymorphism and 
cancer risk was evaluated by calculating the pooled OR and 
95% CI. A Z-test was used to determine the significance 
of the pooled OR, and P<0.05 was considered statistically 
significant. The meta-analysis examined the SNP associations 
in the following models: the co-dominant model (TT vs. CC 
and CT vs. CC); the dominant model (TT+CT vs. CC), and the 



recessive model (TT vs. CC+CT). 

The x 2 -based Q test was used to detect heterogeneity, 
and I 2 index was adopted to measure the extent of 
heterogeneity [41 ' 42 i. A P-value>0.10 [43 ' in the Q test indicates 
a lack of heterogeneity among the selected studies. If no 
heterogeneity was found, the Mantel-Haenszel method was 
used to estimate the pooled OR of all individual studies in 
a fixed-effect model [44] . Otherwise the random-effect model 
(the DerSimonian and Laird method) [44] was used when 
heterogeneity was observed among the selected studies. 

Publication bias was evaluated using the funnel plot 
wherein the standard error of log(OR) for each study was 
plotted against its log (OR). The asymmetry of a funnel plot 
was determined using the Begg's rank correlation test [45] . 
Furthermore, the Egger's linear regression test was used [461 
to measure the association between the mean effect estimate 
and its variance. Subgroup analyses were also performed 
according to the cancer types and ethnicity. If a cancer 
type was evaluated only by one study, this was grouped as 
"other cancer." The ethnicity was categorized into Asian 
and Caucasian descent. Moreover, sensitivity analyses 
were conducted to examine the robustness of the results by 
excluding one study at a time and recalculating the combined 
OR and 95%CI of the remaining studies. 

All of the analyses were conducted using the STATA 
software (version 11.0; StataCorp, College Station, TX) and 
Review Manager (version 5.0.0; The Cochrane Collaboration, 
Oxford, England), with two-sided P-values. 

Results 

Characteristics of the studies 

A total of 39 articles were initially identified through database 
searches using different key words and their combinations. 
After reading the titles and abstracts of the identified articles, 
the studies that did not meet the inclusion criteria were 
excluded. During data extraction, additional articles that did 
not have allelic frequencies available for analysis were also 
disqualified. Eligible studies were retrieved for detailed full- 
text evaluation. Manual searches were also done and one 
relevant article from the references of selected studies was 
identified. Finally, 27 independent case-control studies' 14401 , 
involving 10,435 cases and 12,075 controls, were included 
in the meta-analysis for rsll614913. Figure 1 shows the 
identification and selection of studies. 

The characteristics of the eligible studies are summarized 
in Table 1. As shown in the table, 10 studies involved 
Caucasians and 17 involved Asians. Four studies focused 
on the association of SNPs with lung cancer, 4 with liver 
cancer, 5 with breast cancer, 4 with colorectal cancer, 2 with 
gastric cancer, 2 with esophageal cancer, and the remaining 
6 studies were associated with other types of cancer. Up 
to 12 studies used hospital-based controls and 9 used 
population-based controls. The control sources for the other 
6 studies were not identified clearly. The main genotyping 
method was polymerase chain reaction-restriction fragment 
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Searching in the PubMed, EMBASE, Cochrane library and Web of 
Knowledge database (last update 31 Jan 2012) 





Exclusion of unqualified studies by screening 




titles and abstract (n=13): 


— ► 


studies on non-cancer disease (n=7); 




meta-analysis (m=4); 




Case-only studies (n=2) 



Remaining literatures underwent full tests (n=26) 



+ 



Manual search of the reference lists of retrieved 
studies (n=l) 



Finally included studies (n=27) 



Figure 1. Flow diagram of the identification of studies. 
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Figure 2. miR-196a-2 rsll614913 MAF distribution in different 
ethnicity. 



length polymorphism (PCR-RFLP), which was used in 15 
studies. The other genotyping methods employed included 
MassARRAY, Taqman, DNA sequencing, and high-resolution 
melting analysis. One study was published in 2008, 3 in 2009, 
11 in 2010, and 13 in 2011. 

Quantitative synthesis 

Minor allele frequencies (MAF) were calculated in the 
control groups of each study to estimate the racial difference. 



Figure 2 shows the significant difference (C as the minor 
allele) in MAF distribution of miR-196a-2 rsll614913 T>C 
polymorphism among the controls of the different races: 
Asians and Caucasians. Specifically, the mean MAF in 
17 Asian studies was 0.451±0.013. The mean MAF was 
0.637+0.020 in the 10 Caucasian studies. The distribution of 
rsll614913 genotypes in the control subjects followed the 
HWE in all studies except 2, reported by Mittal et al. [37] and 
by George et al. I 38 ' (Table 1). 

Meta-analysis results 

The associations of the miR-196a-2 polymorphism with the 
risk of developing different cancer types among the different 
racial groups are shown in Table 2. In the meta-analysis of 
all studies, the rsl 1614913 polymorphism was significantly 
associated with the risk of cancer in 4 genetic models (CT 
vs. TT, OR=1.15, 95%CI=1.05-1.27; CC vs. TT, OR=1.23, 
95%CI=1.08-1.39; Dominant model, OR=1.17, 95%CI=1.06- 
1.30; Additive model, OR=1.08, 95%CI=1.01-1.14). 

The tumor type, race, control source, and lab test methods 
were adjusted as confounding factors in the overall analysis. 
Subgroup analysis was performed according to different 
tumor types, racial groups, control sources, and genotyping 
methods. In the subgroup analysis of tumor types, the C allele 
was associated with increased risks for lung cancer, breast 
cancer, and colorectal cancers, but was not associated with 
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those for liver cancer, gastric cancer, 
and esophageal cancers. A significant 
association with lung cancer was 
observed in all genetic models (CC 
vs. TT, OR=1.26, 95%CI=1.07-1.49; 
CT vs. TT, OR=1.15, 95%CI=1.00- 
1.33; Dominant model, OR=1.19, 
95%CT=1.04-1.36; Additive model, 
OR=1.12, 95%CI=1.03-1.22) except 
for the recessive model. Similarly, 
statistically significant associations 
were found for colorectal cancer 
in all models (CC us. TT, OR=1.44, 
95%CI=1.18-1.75; Dominant model, 
OR=1.25, 95%CI=1.06-1.46; Recessive 
model, OR=1.30, 95%CI=1.10- 
1.53; Additive model, OR=1.20, 
95%CI=1.08-1.32). In breast cancer, 
the CC genotype was associated 
significantly with increased risk when 
compared with TT and CT genotype. 

In the subgroup analysis of race, 
statistically significant associations 
with increased cancer risk were 
found among Asians in all genetic 
models. Specifically, the OR (95%CI) 
was 1.27 (1.11-1.46) for CC vs. TT, 1.16 
(1.04-1.28) for CT vs. TT, 1.19 (1.07- 
1.32) for CT+CC us. TT, 1.15 (1.04- 
1.27) for CC vs. CT+TT, and 1.12 (1.05- 
1.20) for C allele vs. T allele. However, 
no significant association was found 
for Caucasians in any of the genetic 
models. 

In the subgroup analysis of the 
control source, the CC genotype 
was associated with cancer when 
the study used a population- 
based control (CC vs. TT, OR=1.24, 
95%CI = 1.11-1.39; Dominant 
model, OR=1.10, 95%CI=1.01-1.20). 
Nevertheless, the CT genotype was 
marginally associated with cancer 
risk in the hospital-based control 
subgroup (CT vs. TT, OR=1.17, 
95%CI=1.00-1.37). 

Subgroup analysis was also 
performed based on different 
genotyping methods. Statistically 
significant associations with increased 
cancer risk were found in the PCR- 
RFLP subgroup in all genetic models 
except for the recessive model. 
Specifically, OR (95%CI) was 1.27 
(1.15-1.41) for CC vs. TT, 1.11 (1.02- 



68 



Zhang et al. Mir-196a-2 Polymorphism and Cancer Susceptibility 



Ye Y (2008) 
Tian T(2009) 
Hu Z (2009) 
Hoffman AE (2009) 
OkuboM (2010) 
DouT(2010) 
Wang K (2010) 
Peng S (2010) 
Qi P(2010) 
Li XD (2010) 
Kim JM (2010) 
Liu Z (2010) 
Srivastava K (2010 
Catucci I (2010) 
Min KT(2011) 
Zhang M (2011) 
Zhan JF(2011) 
Zhou B (2011) 
Hong YS (2011) 
Zhang XW (2011) 
Zhu LJ (2011) 
Chen H (2011) 
Jedlinski DJ (2011) 
Akkiz H (2011) 
Mittal RD (2011) 
George GP (2011) 
Vinci S (2011) 



Meta-analysis estimates, study omitted 
I Lower CI Limit O Estimate I Upper CI Limit 



I- 



1.09 



1.12 



G 
■G- 



1.20 



1.28 



1.32 



Figure 3. Pooled OR of sensitivity analysis. 
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Figure 4. Begg's funnel plot for rsll614913 (CC vs. TT). 



1.21) for CT vs. TT, 1.16 (1.07-1.25) for CT+CC vs. TT, and 1.09 
(1.01-1.18) for C allele vs.T allele. However, no significant 
association was found for genotyping methods other than 
PCR-RFLP in any of the genetic models. 



Test of heterogeneity 

The heterogeneity of the studies was analyzed in the overall 
meta-analysis as well as subgroup analysis. The P values 
in the Q test were also shown in Table 2. The heterogeneity 
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of the study decreased in the subgroup analysis of cancer 
type. If P>0.10, a fixed-effect model was chosen, otherwise a 
random-effect model was used. 

Sensitivity analysis 

In the sensitivity analysis, the meta-analysis was repeated 
by excluding 1 study each time to assess the influence of the 
removed study on the pooled ORs. The corresponding pooled 
ORs were not altered materially for rsll614913 (Figure 3 and 
Table 3). 

Publication bias 

Begg's funnel plot and Egger's test were performed to 
evaluate the publication bias of the study. The funnel plot for 
CC vs. TT of miR-196a-2 rsll614913 polymorphism seemed 
approximately symmetrical (Figure 4). The Egger's test did 
not show any evidence of publication bias (t=0.60, d/=26, 
P=0.553). 

Discussion 

MiRNAs are important post-transcriptional regulators of 
gene expression and are reportedly involved in various 
diseases. Emerging miRNA profiling studies have 
investigated the differences in miRNA expression between 
cancer patients and healthy controls. Various miRNAs 
are upregulated or downregulated in different forms of 
tumors. However, the regulatory effects of miRNA in cancer 
occurrence and development are complicated. Although 
genetic and epigenetic regulation influence miRNA activity, 
the mechanism of miRNAs involvement in carcinogenesis 
remains unclear. Considering the extensive regulation 
of many target genes by miRNA, a single nucleotide 
polymorphism in the miRNAs gene may affect the expression 
and activity of various genes or proteins, which further 
influence the carcinogenic process and tumor progression. 

A previous study found 323 SNPs located in 227 human pre- 
miRNA sequences [101 . Another study integrated 474 human 
miRNA sequences, and found genetic variants located mainly 
in the precursor sequences (10%), whereas rare variants were 
located in the seed region (<1%) [35] . Of these pre-miRNA 
SNPs, the miR-196a-2 rsll614913 T>C polymorphism is the 
most studied SNP This functional polymorphism reportedly 
has an association with susceptibility to various forms of 
tumor including lung cancer, breast cancer, glioma, gastric 
cancer, gallbladder cancer, head and neck cancer, esophageal 
cancer, hepatocellular carcinoma, as well as shortened 
survival time in non-small cell lung cancer. 

The results of the meta-analysis show that people carrying 
the variant genotype C have an increased risk of cancer 
compared to those carrying the wild genotype T. This finding 
coincides with a previous meta-analysis |47 '. The miR-196a-2 
polymorphism may play a role in carcinogenesis. 

Stratification of the analysis in terms of tumor type 
revealed that the miR-196a-2 polymorphism was related 



to lung cancer, colorectal cancer, and breast cancer, but no 
relationship was seen for liver cancer, gastric cancer, or 
esophageal cancer. The differences may be attributed to the 
following: miRNAs may have tissue-specific expression and 
the same miRNA SNP may play different roles in different 
tissues, which leads to different degrees of carcinogenesis; 
the sample size in each subgroup was not large enough to 
have sufficient study power to generate reliable results. 
However, our subgroup analysis showed substantial results 
in accordance with previous studies on various types of 
cancer. Previous meta-analyses has found that rsll614913 is 
associated with colorectal cancer risk 148] , breast cancer risk' 49 ', 
and lung cancer risk [47] , consistent with the results of our 
study. 

In the subgroup analysis by race, an increased risk of 
cancer was found among Asians in all genetic models, but not 
among Caucasians. The distribution of the minor allele in the 
controls of these populations was also evaluated. Significant 
differences in the distribution of C allele were observed. 
Caucasians have a higher frequency of the C allele compared 
with Asians. This difference in allelic frequency may affect 
genetic susceptibility in different racial/ethnic groups. 

In the subgroup analysis in terms of the control source, the 
C allele was associated with cancer susceptibility in studies 
with population-based controls. Moreover, the population- 
based subgroup was less heterogeneous than the hospital- 
based subgroup, which indicated that population-based 
controls are more helpful in reducing the heterogeneity of 
observational studies. 

Finally, the C allele was associated with cancer risk among 
studies that used PCR-RFLP for genotyping. However, in 
subgroups using other genotyping methods, no significant 
association was found. Different genotyping methods may 
have different capacities for detecting gene polymorphism, 
which may influence the results. The laboratory test method 
may be one of the confounding factors in the overall analysis. 
After adjusting for genotyping method, the rsll614913 C 
genotype was associated with cancer risk. 

Currently, two possibilities explain the association of 
the miR-196a-2 polymorphism with cancer risk. First, the 
rsll614913 variant may increase mature miR-196a expression 
and enhance target mRNA binding |13 L Zhan et al. [26] reported 
that the miR-196a expression level among C allele carriers 
was higher than that in those carrying the TT genotype in 
colorectal cancer. Similarly, Li et al. found increased miR-196a 
expression among hepatocellular carcinoma patients with 
Hepatitis B virus infection who carried the C allele [22] . This 
SNP is associated with G 2 cell cycle delay, which plays a vital 
role in carcinogenesis' 171 . 

The current meta-analysis has several limitations. First, 
the cancer types included in the analysis were limited. 
Although the polymorphism and cancer risk in all reported 
studies were analysed, these studies mainly focused on 
several tumor types, which may affect the risk estimation. 
Second, some of the study controls were not from healthy 
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populations. Selection bias could have occurred, which may 
have confounded the results. Third, the lack of individual 
information inhibited the calculation needed to adjust the 
ORs and analyze the gene-environment interaction. 

In conclusion, miR-196a-2 polymorphism is associated 
with cancer risk, especially lung cancer, colorectal cancer, 
and breast cancer. Asians carrying the miR-196a-2 variant 
genotype are more susceptible to cancer compared with 
Caucasians. The SNP in miR-196a-2 may be a key factor in 
carcinogenesis. Furthermore, well-designed studies involving 
various ethnic populations and more cancer types are needed 
to confirm the results. Functional studies are needed to clarify 
the mechanisms that show the effect of miRNA SNPs on 
cancer development. 
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